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ABSTRACT 

Motivated by recent observations that show increasing fractional linear polariza- 
tion with increasing wavelength in a small number of optically thin jet features, i.e. 
"inverse depolarization", we present a physical model that can explain this effect and 
may provide a new and complementary probe of the low energy particle population 
and possible helical magnetic fields in extragalactic radio jets. In our model, structural 
inhomogeneities in the jet magnetic field create cancellation of polarization along the 
line of sight. Internal Faraday rotation, which increases like wavelength squared, acts to 
align the polarization from the far and near sides of the jet, leading to increased polar- 
ization at longer wavelengths. Structural inhomogeneities of the right type are naturally 
produced in helical magnetic fields and will also appear in randomly tangled magnetic 
fields. We explore both alternatives and find that, for random fields, the length scale 
for tangling cannot be too small a fraction of the jet diameter and still be consistent 
with the relatively high levels of fractional polarization observed in these features. We 
also find that helical magnetic fields naturally produce transverse structure for inverse 
depolarization which may be observable even in partially resolved jets. 

Subject headings: galaxies: active — galaxies: jets — radiation mechanisms: non- 
thermal — radio continuum: galaxies — quasars: general — BL Lacertae objects: 
general 



1. Introduction 



Faraday rotation, the rotation of the plane of linearly polarized radiation propagating through 
a magnetized plasma, has long been thought not to play an important role internal to extragalactic 
radio jets, in large part because the depolarization and diverge nce from a pure A 2 rotation law 
expected due to large amounts of internal rotation (|Burnlll966T) have not been cl early observed; 
however, detecting this behavior in real sources can be difficult (jCioffi Jonesll 19801 ). Observations 
of circularly polarized radiation have suggested that Faraday rotation may indeed operate internal 
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to jets to drive th e conversion of linear polarization into circular (e.g. IWardle fc Homanl 120031 ; 
Homan et al.ll2009l ). Additionally, observations of variations in Faraday rotation and transverse 
structure, or gradients, in rotation measure have indicated that some Faraday screens lie very close 
to the jet, perhaps in a boundary layer around it that may contain helical magnetic fields (e.g. 



2011 



Asa da et al.l 120021 ; iGabuzda et al J 12004 : IZavala & Tavlorll2005l ; iGomez et al.ll201ll ; iKronberg et al 



). The observation and interp retation of transverse rotation measure gradients as evidence for 
helical fields remains controversial (j Taylor &; Zavalall 20101 : iBroderick &; McKinneyil2010l ) . Given the 
central role helical fields are believed to play in collimating a nd accelerating je ts near their origin 



from the super massive black hole-accretion disk system (e.g. iMeier et al.ll200lf), it is important to 



pursue additional observational signatures of he l ical magnetic field s in jets (jLvutikov et al 



Marscher et al. 



2008 



Clausen-Brown et al 



2011 



Porth et al 



2011 



2005 



An intriguing phenomenon has recently been observed in multiple optically th in jet features 



i n the blazars 3C 273, 3C 454.3, and a single jet feature in the blazar 1514—241. iHovatta et al 
(|2012. submitted! ) have measured increasing polarization with increasing wavelength in their Very 
Long Baseline Array observations of these jet features at AA 2.0, 2.5, 3.6, and 3.7 cm. The observed 
i ncrea ses range from factors of ~ 1.5 up to ~ 3 times over this wavelength range. iHoman et al 



(|2002l ) saw a similar increase in fractional polarization of up to a factor of two in a jet feature 
in 3C 120 between AA 1.3 and 2.0 cm. The jet features involved in these cases are relatively 
strongly polarized at their maximum wavelength, with fractional polarization ranging from 5 to 
more than 20 percent. The effect can apparently vary with epochs separated by a few months, 
although it is clearly observed for two epochs spanning three months in two features in 3C 273 



(Hovatta et al 



2012. 



3C120 ( Homan et al 



200 



3T 



sub mitted! ) and over four epochs spanning six months in the jet feature in 



The sense of this effect is the opposite of that normally expected for Faraday depolarization 
where fractional polarization decreases with increasing wavelength, and we therefore term it "inverse 
depolarization". Increasing fractional polarization with wavelength might naturally arise in the 
unresolved, optic ally thick cor e region where one is observing different parts of the jet at different 



wavelengths re.g. lKonigllll98lh . but the jet features showing this effect are (1) well separated from 



the core in the optically thin part of th e jet and (2) reasonably isolated from other strong jet 
features (IHovatta et al.l 12012. submitted! ). In this letter we propose an alternative explanation 
which combines structural inhomogeneities in the jet magnetic field along the line of sight with 
internal Faraday rotation to produce a inverse depolarization effect in optically thin jet featured- 
If our explanation is correct, observed inverse depolarization is evidence of internal Faraday rotation 
in jets and may prove to be a powerful and complementary tool to rotation measure studies of jet 
magnetic fields and particles, including helical field geometries. 



*A similar effect due to Faraday conversion of linear to circular polarization was seen in simulations by 
Baldwin fc Homan ( 20091 ) at modest optical depths, r ~ 0.1 — 1.0, while the mechanism proposed here can operate 



in optically thin regions and requires only internal Faraday rotation. 
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Our model for inverse depolarization is presented in §2 and applied to both helical field geome- 
tries and randomly tangled magnetic fields. §3 discusses these results and the prospects for using 
this effect to study the magnetic field geometry of extragalactic jets. 



2. Model for Inverse Depolarization 



In the subsections that follow, we present a simple model for generating inverse depolarization, 
i.e. increasing polarization with increasing wavelength, in radio jets at low optical depth, t<1. 
Our model combines internal Faraday rotation with structural differences in the magnetic field 
between the far and near sides of the jet. If the magnetic field structure naturally leads to low net 
polarization at short wavelengths, internal Faraday rotation can act to align polarization from the 
far side of the jet with polarization produced at the near side of the jet, reducing the cancellation 
between them and leading to increased net polarization at longer wavelengths. 

We characterize the amount of internal rotation at a given wavelength by the Faraday depth, r/, 
which is proportional to A 2 . Internal to the jet, Faraday depth may be due to either a population of 
thermal electrons or the low-energ y end of the relativistic particle distribution (e.g. I Jones &: O'Dell 



1977 



Huang fc Shcherbakov 



2011 



and we discuss these contributions in N2.5 



2.1. Linearly Polarized Transfer at Low Optical Depth 



The radiative transfer equations for Stokes Q and U at low optical depth a re given by the 
following express ions, assuming the projected magnetic field is purely east-west (| Jones fc O'Dell 
19771 : IJonedll988l ). Note that in the low optical depth limit, we have ignored the coupling of Stokes 
I and V to Q and U 



The solutions to these equation^, including an initial polarization, Pq at angle xo are 



(1) 



Q 



eJ 



U 



{rf/rf + 1 
eJ 

(t//t) 2 + 1 



1 — e T cos Tf + ( — J e r sinTy + P$e T cos(2xo + Tf) 



— — e T cosTy) — e r sin7y + P^e T sin(2xo + Tf) 



(2) 



2 The solutions to the general form were obtained with WolframAlpha, h ttp:/ /www.wo lframalpha.com/ and 
checked against numerical integration of the full equations of radiative transfer ( Homan et al. 20091 ) . 
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where r is the optical depth and Tt = Cy T is the Faraday depth and is proportional to A 2 (e.g. 

Jones & 0'Del3ll977l : Icioffi fe Joneslll98oh . Note that at low optical depth, Stokes I = Jr. The 
final Stokes Q and U are used to compute the fractional polarization, m 
electric vector position angle (EVPA), \ = 0.5 



^Q 2 + U 2 /I, and 
tan -1 ([//(3) of the emergent polarization. 



2.2. Inverse Depolarization in a Simple Two-Cell System 

We apply these expressions to a "two-cell" system along the line of sight. Each cell has an 
internal Faraday rotation of depth rt /2 for a total Faraday depth of Tf for the system. In addition 
to their common line-of-sight magnetic field, both cells also have a uniform magnetic field, B±, in 
the plane of the sky. The field in the two cells differs only by a rotation on the sky of A<p. The B± 
in the cell closest to the observer is taken to be purely East- West, while the field in the furthest 
cell differs in orientation by an offset angle, A<p, Northward from East. Note that positive Faraday 
depth will tend to rotate the x of the emitted radiation tan North to Eaafl 

Polarized radiation generated in the furthest cell must travel through the nearest cell to reach 
the observer. In the absence of internal Faraday rotation, the offset angle A<p will produce some 
cancellation between the polarization emitted from the far and near cells. For A<p = 90°, this 
structural cancellation would be complete at short wavelengths. 

Figure 1 shows the emergent polarization as a function of tj (oc A 2 ) for a range of A(j) offsets. 
For A<f> = 0, 1 Burn! (119661 ) depolarization in a slab geometry is reproduced. At larger offsets, partial 
cancellation of the emergent polarization is clear at Tf = 0, and for < A(j) < 90° the internal 
rotation initially begins to reverse this structural depolarization up to a certain limit until 7/ 
becomes too large. For A<p > 90° the internal rotation initially increases the depolarization as tj 
increases but later can also go into a region of 'inverse depolarization'. 

Figure 1(b) shows that for simple field geometries, such as the two cell system shown here, the 
observed EVPA for internal rotation is linear in A 2 with y = Xo + r //4 with additional 90° flips 
when the polarized flux passes through zercQ (e.g. Burn 1966 ). More complicated geometries which 
include a range of tj value s along diff erent lines of sight are expected to diverge from a A 2 law 
after total rotations ^ 45° (|Burnlll966l ). although the total rotation at which the non-A 2 behavior 
becomes apparent will vary from case to case (see the lower panels of Figure 2 as examples). 



Negative Faraday depth would simply result in a rotation of \ in the opposite direction. 

4 Note that for purely external rotation in a uniform screen t/ is also proportional to A 2 , but x ~ Xo + T f /2. 
However, that is not the case we consider in this paper. 
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2.3. Inverse Depolarization in Helical Fields 

To evaluate the efficiency of this mechanis m in a purely helical held geometry , we numerically 



solved the full equations of radiative transfer ([Jones &: O'Delll 119771 : |Joned ll988) for a variety of 



combinations of toroidal (ft) and uniform (f u ) field components in an optically thin cylinder. The 
uniform field is poloidally directed along the jet axis and has constant magnitude throughout the 
cylindrical cross-section of the jet. The magnitude of the toroidal field component grows linearly 
with distance from the jet axis from to ft at the outside edge of the cylinder, consistent with field 
generated by a current which is carried uniformly in the jet volume. A mor e detailed descriptio n 



of the radiative transfer simulation and magnetic field parameters is given by lHoman et al.l ([2009 ). 
The simulated jets are relativistic, with bulk Lorentz factor V = 10, and viewed at one-half the 
optimum angle for superluminal motion 9 = 1/2F This intermed iate viewing angle is the most 



probable in a flux-density limited sample (jLister &: Marscherlll997l ); however, we note that a wide 



range of viewing angles is possible and may strengthen or weaken the inverse depolarization effect 
in various cases. Internal Faraday depth is provided by the low end of the relativistic particle 
spectrum^!, resulting in (r/) ~ 2A 2 , averaged over all lines of sight, for the simulations shown in 
Figure 2. 

Figure 2 plots the integrated polarization from the entire jet cross-section as a function of A 2 for 
a variety of combinations of f u and ft- Even in these integrated results, the inverse depolarization 
effect is clearly seen, for at least some range of wavelengths, in half of the cases shown. Note that 
there is not necessarily any correlation in the integrated results between the amount of inverse 
depolarization and the observed Faraday rotation (given in the lower panel of Figure 2). This lack 
of correlation in the integrated results is primarily due to variation in both Faraday depth and field 
order across the jet. 

In a helical magnetic field, the offset angle, A(j>, between the back and front of the jet varies 
with location along the cylindrical cross-section. Additionally, the internal Faraday depth, 7y, is 
produced by the line-of-sight components of both the poloidal (uniform) and toroidal fields. The 
contribution to ry from the toroidal field will also vary with location along the cylindrical cross- 
section, changing the sign of its contribution at the mid-point of the jet. The details of how 
these factors combine depend on a variety of factors, but the important point is that in a helical 
magnetic field, the inverse depolarization effect can vary strongly as a function of location along the 
jet cross-section. Detecting and modeling this variation requires sensitive observations at several 
wavelengths; however, in Figure 3 we demonstrate that the effects could be observable even with 
only two beam- widths of resolutions across the jet. 

Figure 3 shows a partially resolved, cylindrical jet cross-section in relative total intensity (solid 
curve) and fractional polarization at a range of wavelengths (symbols). Even with only two beam- 



We note that the same results would be obtained by introducing thermal particles to do the rotation, as it is 
only the amount of rotation along each line of sight that determines the emergent polarization. 
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widths of resolution, the inverse depolarization effect clearly varies with location along the jet, being 
extremely strong at the bottom of the jet (left of figure) and turning to normal depolarization at 
the top of the jet (right of figure). This is a particular case, and, as noted above, the details will 
vary depending on the combination of poloidal and toroidal fields, jet viewing angle, and average 
Faraday depth. It is interesting to note that a jet viewed perpendicular to its axis in the fluid frame, 
9 = 1/r, can still have a significant inverse depolarization effect, even though it has no component 
of poloidal field along the line of sight. This can happen if the field crossing angle is A0 ~ ±90° 
above and below the jet axis. Although the ry provided by the toroidal field reverses across the 
axis, either sign of Tf will increase the net polarization from near zero giving inverse depolarization 
on both the top and bottom of the jet. 



2.4. Inverse Depolarization in Random Fields 

The structural inhomogeneities that are necessary for inverse depolarization can also be pro- 
duced stochastically in a randomly tangled magnetic field. The effect should be strongest in those 
cases where structural cancellation leads to low net polarization so that internal rotation can ef- 
fectively increase the fraction of polarization with wavelength. We have simulated this scenario 
with a N x N x N cube of cells, each with a randomly oriented magnetic field. The results for 
N = 5 and N = 15 are plotted in Figure 4 where we show the percentage of cases with significant 
inverse depolarization between two wavelengths (top panel) and the average polarization of those 
cases showing the inverse depolarization (bottom panel). Significant inverse depolarization occurs 
for 20 — 30% of cases for both N = 5 and N = 15, with and without a modest "shock" applied to 
the entire cube (see caption), for a range of mean Faraday depths. However, the mean fractional 
polarization in the N = 15 cases with inverse depolarization is very low, ~ 1 — 2%. The N = 5 cases 
are larger, ~ 5 — 8% polarization, and more consistent with the levels of polarization seen in t he 



jets showing this effect, although levels up to 20% are observed (jHovatta et al.l 12012. submitted! ) . 



2.5. Internal Faraday Depth 



Internal to the jet, Faraday depth may be due to either a population of "cold" thermal electrons 
with number density, n c , or the low-energy end of the relativistic particle distribution with number 
density, n r . The relative contributi ons to the Faraday depth from the cold and relativistic particles 
are given by I Jones & O'Delll (|1977l ): 



T f 
r (c) 



2a 



a + 3/2 ln 7i 



71, 



n, 



(3) 



a + 1 Yi 

where a is the optically thin spectral index defined by S oc v~ a . The relativistic contribution as- 



-7- 



sumes a power-law distribution of particles with dn r = Kj ( 2a+1 )(f7 with a lo wer cutoff ^j. For cold 



parti cles alone, we have the following familiar expression (in cgs units) (e.g. iRvbicki &; Lightman 



197a ). 



rf = & \ 4 / n c B • dl (4) 

J 7rm z c 4 (1 + z) z J 

Note that A in the above expression is the observed wavelength. Internal to the jet, the rotating 
plasma sees a much longer wavelength because the jet emission is Doppler boosted in the observer's 
frame by factor, 5/ (1 + z). It is interesting to note that, for typical blazars with 8 ;> 10, low-energy 
particles internal to the jet (whether cold or relativistic) can be hundreds of times more effective 
at Faraday rotation than rotating plasma external to the jet which is not moving relativistically in 
relation to the observer (assuming equivalent field strengths and path lengths). 

The inverse depolarization effect proposed here requires only modest Faraday depths of order 
Tf ~ a few, and indeed large Faraday depths would effectively depolarize jets, which is not ob- 
served in blazars. Suppression of large internal Faraday depths requires some combination of the 
following conditions: (1) low number density of the rotating particles, both thermal and low-energy 
relativistic, (2) pairing of low-energy relativistic electrons with low-energy relativistic positrons 
without a significant thermal component, or (3) a large number of magnetic field reversals along 
the line of sight. Study of this effect which implies modest but not large internal Faraday rotation, 
and provides constraints on the magnetic field order, may yield new information on the low-energy 
particle population of relativistic jets. 



3. Summary and Conclusions 



We have shown that observation s of inverse depolarization in optically thin jet features (jHovatta et al. 
2012. submitted ; Homan et al.ll2002 ) , can potentially be explained by a simple physical model which 
combines structural inhomogeneities in the magnetic field along the line of sight with internal 
Faraday rotation. Structural inhomogeneities of the right type are naturally produced by helical 
magnetic fields. Random magnetic fields can also produce this effect, but must be tangled on long 
length scales in the jet to be consistent with the modest to high levels of fractional polarization 
observed in these features. R eal jets likely contain some com bination of randomly tangled and 



ordered field components (e.g. iHughes 



2005 



Laing et al.ll2006l ); however, it is intriguing that 3 of 



the 4 jets which appear to show this effect also have transverse rotation measure gradients, a po- 



tential signature of helical magnetic fields (jAsada et al. 



2008 



2002 



Zavala Sz Tavlorl 120051 ; iGomez et al 



Hovatta et al.l 12012. submitted! ), although for 3C 120 the evidenc e favors an interact ion with 



the external medium as the source of the rotation measure gradient (jGomez et al.l 120081 ). In the 
case of inverse depolarization produced by helical magnetic fields, we have demonstrated that the 
transverse structure of inverse depolarization could be measured even in partially resolved jets, 
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and the recently expanded bandwidth of the Very Long Baseline Array may allow highly sensitive, 
multi-wavelength observations of transverse jet profiles to test these predictions in the near future. 

The author would like to thank John Wardle, Talvikki Hovatta, Matt Lister, Margo Aller 
and the other members of the MOJAVE collaboration for helpful discussions. This work has been 
supported by NSF grant AST-0707693 
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T f (cc X 2 ) 

Fig. 1. — Plot of fractional linear polarization (top panel) and EVPA (bottom panel) as a function 
of Faraday depth, t/ oc A 2 , for a simple two-cell system along the line of sight. Each line shows the 
emergent polarization for a different offset angle, A<f>, for the projected magnetic field between the 
far and near cells. 
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Fig. 2. — Integrated polarization as a function of A 2 for cylindrical jets with pure helical magnetic 
fields. Jets with various combinations of uniform magnetic field, along the jet axis, f u , and toroidal 
magnetic field, ft are plotted. All of the jets here are identical except for their (f u , f t ) combination, 
and all are at an angle 6 = 1/2T to the observers line of sight, inside the optimum angle for 
superluminal motion. 



,x =4.0 




Jet Cross Section Location (beams) 



Fig. 3. — Fractional polarization, m, as a function of transverse jet position for a partially resolved 
jet where the jet diameter is twice the FWHM of the observing beam. The (f u , ft) = (0.4,-0.6) 
case from figure [2] is plotted for several values of A 2 . The single solid line represents the Stokes / 
intensity normalized to a maximum of 1.0 . Note that the asymmetric profile in Stokes I is similar 
to that predicted by IClausen-Brown et al.l (|201ll ) for a helical magnetic field viewed at 9 = 1/2T. 
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Fig. 4. — Top panel shows the percentage of 1000 runs where the fractional polarization ratio, ffiratio 
between A 2 = 2.0 and A 2 = 1.0 is > 1.5 indicating significant inverse depolarization. Integrated 
results from cubes of random field of dimensions 5x5x5 and 15 x 15 x 15 are plotted with and 
without a moderate shock where unit length has been shortened to length k along the jet axis 



(jWardle et al.lll994l ). Each case is plotted for a variety of mean Faraday depths, < tj >. The 



bottom panel shows the mean fractional polarization at the longest wavelength for the runs where 

Horatio ^ 1-5. 



